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In a recent paper we derived an expression for the replicated free energy of a liquid of hard 
spheres based on the HNC free energy functional. An approximate equation of state for the glass 
and an estimate of the random close packing density were obtained in d — 3. Here we show that the 
HNC approximation is not needed: the same expression can be obtained from the full diagrammatic 
expansion of the replicated free energy. Then, we consider the asymptotics of this expression when 
the space dimension d is very large. In this limit, the entropy of the hard sphere liquid has been 
computed exactly. Using this solution, we derive asymptotic expressions for the glass transition 
density and for the random close packing density for hard spheres in large space dimension. 



I. INTRODUCTION 



The study of amorphous packings of hard spheres is relevant for a large class of proble ms, including liquids, glasses, 
colloidal dispersions, granular matter, powders, porous media, etc. 0,UM Lm^ta&y,UilUJ,L2, 131. Nevertheless, 
the question whether a glass transition exists for a liquid of identical hard spheres in finite dimension is still open. 

Recently, a quantitative description of the glass transition in structural glasses has been obtained by means of 
the replica trick El El El El Il8l Il9| . The latter is applied in the context of a variational approximation - the 
hypernetted chain (HNC) approximation - that leads to a suitable free energy functional for simple liquids [2I . 
This method was successfully applied to Lennard- Jones systems |l4l Il5l lltil Il7| and, more recently, to hard spheres 
in space dimension d = 3 El E3 ■ Quantitative estimates of the glass transition temperature (or density) and of the 
equation of state of the glass have been obtained in this way. Moreover, for hard spheres an estimate of the random 
close packing density, i.e. the maximum density of the amorphous configurations of the system, has been obtained as 
the value of the density where the pressure of the glass diverges EH • 

In this approximation the glass transition turns out to be similar to the 1-step replica symmetry breaking (1RSB) 
transition that happens in a class of mean-field s pinglas s models and indeed the replica strategy described above was 
inspired by the exact solution of these models [2H2! |2J, |2f| . However, for finite dimensional models with short range 
interactions the picture emerging from the replica-HNC approach should be modified by non-perturbative activated 
processes (for a detailed discussion see e.g. [H|2l|2l|). The results obtained in d = 3 have often been considered as 
a kind of "mean-field" approximation. 

Activated processes should become less relevant on increasing the space dimension. It is then natural to ask whether 
this approximation describes better and better the true properties of the system (eventually becoming exact) when 
the space dimension is large as it is well known that the mean field approximation becomes exact in the d — > 00 
limit. Moreover, the study of sphere packings in large space dimension is relevant for information theory and lot of 
effort has been devoted to finding the densest packing for d — ► 00 [2^. l30l|. Despite this effort, only some not very 
restrictive bounds have been obtained, and it is still unclear whether the densest packings for d — > 00 are amorphous 
or crystalline. 

In this paper we will: i) show that in d — > 00 the replica approach predicts the existence of a glass transition density 
pK and compute the value of this density; ii) compute the maximum density p c of the amorphous packings (or glass 
states) in d — > 00. Unfortunately the value of p c we find is within the current bounds for crystalline packings so we 
cannot address the problem whether the densest packings are amorphous or not. 

The paper is organized as follows: first we briefly outline the method, without entering in details as they are 
discussed in EH; then we present an improvement of the theory of EH) i- e - we show how to construct the small cage 
expansion without using the HNC approximation; finally we discuss the limit d — > 00 of the theory. 
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II. THE METHOD 



In the replica method one considers a liquid made of m copies of the original system, with the constraint that each 
atom of a given replica must be close to an atom of the other m — 1 replicas, i.e. that the replicated liquid must 
be made of molecules of m atoms, each one belonging to a different replica. It was shown that this trick allows to 
compute all the properties of the glass, including the size of the cag es, the vibrational entropy, etc., provided one is 
able to perform the analytical continuation to real m < 1, see |l4 ITU Il6l l25| for a detailed discussion. Thus the 
problem is to compute the free energy of the replicated system for integer m and to continue the expression for real 
m. 

This strategy has been first tested on mean field spin glass models |2^ , and then applied to systems of particles 
interacting through a Lennard- Jones like potential 14, 16, 17] . To compute the replicated free energy for a system of 
particles, the idea was to start from the standard HNC free energy for a molecular liquid [2lJ and expand it in a power 
series of the "cage radius", that represents the amplitude of the vibrations of the particles in the glass state 15]. The 
method was successful but could not be extended straightforwardly to hard spheres, because at some stage it was 
assumed that vibrations were harmonic, an approximation that clearly breaks down for hard core potentials. 

The small cage expansion of the replicated HNC free energy was worked out in ^^] . It was shown that the theory is 
in reasonable quantitative agreement with numerical data even if the HNC approximation is a very poor description 
of the hard spheres liquid. In next section we will show that the result of |l9j for the replicated free energy, obtained 
starting from the HNC free energy functional, can be derived from the full diagrammatic expansion, without the need 
of neglecting any class of diagrams. 



III. SMALL CAGE EXPANSION 



Our approach will be based on the standard diagrammatic approach (virial expansion). It seems to us that similar 
results could also be obtained by a direct computation, however we think that it is more instructive to use the more 
familiar diagrammatic approach. 

We start from the diagrammatic expansion of the canonical free energy as a function of the single-molecule density 
p(x) and of the interaction potential between molecules. It can be obtained from the grancanonical partition function 
of the system via a Legendre transform I2j| . Calling x = (xi, . . . ,x m ) the coordinate of a molecule, the expression 
for the replicated free energy functional is |20ll2l| (setting from now on f3 = 1): 



*[>(x),/( X) y)]= |rfxp(x)[logp(x)-l]- * V^ + IZj + IZl 



(1) 



where the black circles represent a function p(x) and the lines represent a function /(x, y) = e~ ip< ' x ' y ' > — 1, and 
</?(x, y) = ^2 a (p{\x a — Va\) is the interaction potential between two molecules, with ip{r) — oo for r < D and 
otherwise, the usual hard sphere interaction. The function /(x, y) is equal to —1 if \x a — y a \ < D for at least one pair 
a of spheres, and vanishes otherwise. The information on the inter-replica interaction that is used to build molecules 
is encoded in the function p(x). 

The generic diagram of the series expansion of 4" represents an integral of the form 



dxidx2dx3p(xi)p(x2) / o(x 3 )/(xi, x 2 )/(x 2 , x 3 )/(x 3 , xi) 



(2) 



where S is the symmetry number of the diagram [20 
For the one-molecule density we make the ansatz 
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p(x) = p \ dX\\ - — = p / dXp(u) , p(u) = TT ; 3 (3) 

where p = N/V is the number density of molecules, X is the center of mass of the molecule and u a = x a — X . This 
allows to compute exactly the ideal gas term of the free energy: 

1 f d d 

— j dxp(x)pogp(x)-l] = logp-l + -(l-m)log(27rA) + -(l-m-logm) . (4) 



We want to perform a power series expansion in \[A of the interaction term. 



3 



At the zeroth order, the functions p(u) are products of delta functions, p(u) = Y[ a H u a), so that all the coordinates 
x a are equal to X. This gives a contribution 



Zl\ = |- J dX 1 dX 2 dX 3 F(X 1 ,X 2 )F(X 2 ,X 3 )F(X 3 ,X 1 ) 



(5) 



where F(X,Y) = —9(D — \X — Y\) is represented by a dashed line and a vertex i without a black dot represents 
the constant p and an integration over the corresponding Xj. This is exactly the contribution of the diagram to the 
usual non-replicated free energy. Thus, at the zeroth order the molecular free energy simply reduces to the free energy 
T\p, F(r)] of a liquid made by the center of masses of the molecules, with density p and an hard core interaction. 
To discuss the first order correction in \J~A let us first insert Eq. into J2J): 



A 



j ' dX l dX 2 dX 3 du 1 d\i2du 3 p{\i l )p(u 2 )p{\\ 3 )f{X 1 +Ui,X2+\i2) , (6) 



S 

using the notation X + u = (X + u\, ■ ■ . ,X + u m ) and 

f(X 1 + u 1 ,X 2 + u 2 ) = e -S.KS-x 2 +« 1 .-t ls ») _ 1 ( 7 ) 

For small A the u are small too, u ~ \f~A. The function / above has the property that, if \X\ — X 2 \ differs from D by 
a quantity yC4, it is independent of Ui, 112: in fact it is a constant equal to —1 if \X\ — X2I < D and otherwise. 
That is, if \Xi — X 2 \ is not close to D, one has f(Xi + Ui, X 2 + u 2 ) = F(Xi, X 2 ) and the diagram does not give any 
contribution apart from the zeroth order discussed above (recall that J rfup(u) = 1). 

Thus, the corrections in \/~A due to a given diagram come from the region of the integration space where at least 
two of the coordinates Xi, Xj connected by a link have distance D, \Xi — Xj\ ~ D + 0(VA). Let us call such a pair 
a "singular pair" , and the link connecting them as a "singular link" . The regions of the integration space of the Xi 
where n pairs Xi, Xj are singular have a total volume ~ (\/A) n . Thus the only contribution at order yG4 comes from 
n = 1, i.e. only one singular pair. 

In all the non-singular links we can replace /pQ + u i; Xj + Uj) — > F(Xi,Xj) and we get 

A-U 

dxidx2dx3p(xi)p(x2)p(x 3 )/(xi, x 2 )/(x 2 , x 3 )/(x 3 , xi) 

J du 1 du 2 p(u 1 )p(u 2 )/(X 1 + Ul ,X 2 + u 2 ) - F(X l7 X 2 ) 



/--\ + I dX 1 dX 2 dX 3 



S 



F(X 2 ,X 3 )F(X 1 ,X 3 ) 
(8) 



where the wiggly line represents the function 

Q(X U X 2 ) = J du 1 du2p(u 1 )p(u2)f(X 1 +u u X2 + U2)-F(X 1 ,X 2 ) , (9) 



which is different from only if \Xi — X 2 \ ~ D + 0(v A). Thus for each diagram the correction is obtained by replacing 
one F-\mk with one Q-link in all the possible non-equivalent ways. This is equivalent to taking the derivative of the 
zeroth order diagram with respect to F(Xi, X 2 ), multiplying by Q(X±, X 2 ) and integrating over X±, X2. Summing 
the contribution of all the diagrams we then obtain the first order correction to the replicated free energy: 

A*[p(x),/(x,y)] = J dXidX 3 — Q(X 1 ,X 2 ) . (10) 

where the derivative with respect to F(X\, X 2 ) of the free energy is taken at constant p. Recalli ng t hat the canonical 
free energy is the Legendre transform of the grancanonical free energy, it easy to prove that |2(J, |2l| 

8T SlogZac =-£y(x 1 ,x 2 ), (11) 



SF{Xx,X 2 ) SF(X 1 ,X 2 ) 2 



where the function Y(r) = e v ^G(r) is continuous for hard spheres. 
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The calculation of the function Q(X\ , X2) = Q (\X\ — X2 | ) was already done in and is reported in Appendix [X] 



Q(r) = 2VAQ m -£ d (D)5(r - D) , 



(12) 

d/2jjd-l 



mm 



where Q m is an analytic function of m, with Q m ~ 0.638(1 — to)+o((1 — to) 2 ), see Appendix 1X1 and Yj d {D) = — 
is the surface of a d-dimensional hypersphere of radius D. 

Substituting Eq.s (JTTJ, (jT^jl in Eq. (JTOJ, and adding the correction coming from the ideal gas term, Eq. (@J, we 
obtain 



*[p(x),/(x,y)]=JF(p)+iV 

= Hp) + n 



~ (1 — m) log(27rA) + ~ (1 — m — log to) 
I (1 - to) log(27rA) + ~ (1 - to - log m) 



/•OO 

p 2 VAQ m E d (D)V / drS(r - D)Y(r) 
Jo 

NpVAQMD)Y(D) , 



(13) 



where T(p) is the free energy of the non-replicated liquid and Y(D) = G(D) is the value of the pair correlation 
function at contact. This is the same result obtained in from the expansion of the HNC free energy, but here it 
was obtained without need of the HNC approximation. This is important because the HNC approximation gives very 
poor results for hard spheres systems. 

Recalling that T{p) = —NS{p), where S(p) is the equilibrium entropy per particle of the liquid, we obtain the 
replicated free energy by optimizing with respect to A in Eq. (|13l) : 



$(to, p) = i min *[p(x), /(x, y)] = -S(p) + -(1 - to) log[27rA*] + - (to - 1 - log to) 



y/A*(m) 



where Vd(D) = 2 J r ^ d ^ is the volume of a sphere of radius D. The function A*(m) has the meaning of a "cage radius" 
as discussed above. This result holds in any space dimension d. 

In d = 3 we used the Carnahan-Starling expression pfjj for the entropy S(p), which reproduces very well the 
numerical data for the equation of state of the hard sphere liquid. This was done in |l9j on a phenomenological 
ground, but is fully justified by the present derivation in which no reference to the HNC approximation is done. It 
was shown that Eq. I|14|) predicts a glass transition density px which is in good agreement with numerical results. 
Moreover, in the glass phase the cage radius decreases with the density and reaches at a value p c which is then the 
maximum allowed density for an amorphous state, i.e. the random close packing density. For p — > p c the pressure of 
the glass diverges. It was also shown that the average number of neighbors of a given sphere at p c is equal to z = 2d. 
A detailed discussion and a comparison with numerical data can be found in |19| . Here we will discuss the predictions 
of Eq. (|T1J| in the limit d — > 00. 
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D 



(14) 



Q m pV d (D)Y(D) 



IV. ENTROPY OF THE LIQUID IN THE LIMIT OF LARGE DIMENSION 

The problem of computing the entropy of the hard sphere liquid for d — > 00 was addressed in [3lL l32| , where the 
same result was obtained in two independent ways. In [3ll | it was shown that the ring diagrams dominate the virial 
series order by order in p for large d, and the entropy was computed by a resummation of these diagrams; in [32| 
simple equations for the pair correlation function g(r) were introduced, and solved in the limit d — > 00. In both 
cases it was found that S(p) is given by the ideal gas term plus the first virial correction (i.e. by the Van der Waals 
equation). 

The equations introduced in |32| are indeed the minimal requirements for a pair distribution function g(r). If D is 
the sphere diameter, p = N/V is the density, h(r) — g(r) — 1 and h(q) is its Fourier transform, one has: 

'g(r) > , 

< g(r) = for r > D , (15) 
S(q) = 1 + ph(q) > , 

The first condition comes from the fact that g(r) is a probability (the probability of finding a particle at distance r given 
that there is a particle in the origin), the second from the fact that two particles cannot be at distance smaller than 
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D (due to the hard core repulsion). The third condition is easily obtained by proving that S(q) = 1 + ph(q) = (\p q \ 2 ), 
where p q is a Fourier component of the density fluctuations |20| . 

In [23| a solution of E o.s for d — > 00 was found. We set the sphere diameter D = 1 (defining Vd = V^(l) and 
Y = Y(l)) and following [32j we define d = 2N + 3, the reduced density p = pVd and 

_logp_21ogp 

Pl = ^r = i^j- (16) 

Note that p is related to the packing fraction <f> = pVd(l/2), i.e. to the fraction of volume covered by the spheres, by 
<\> = 2~ d ~p. It is found that for p < 1 the solution of Eq.s (|15|l is simply g(r) = 0(r — 1), while for p > 1 it has the form 
g(r) = 9(r — 1){1 + exp[— Nh\(r)]} with h\{r) > 0, i.e. it is given by the step function plus an exponentially small 
correction which can be explicitly computed, see |32j. The pressure is then found to be: 




for p < 1 , ( 7) 



where 



and p c the solution of 



K(p) - log(l + Vl-P 2 ) - Vl-P 2 - logp , (18) 



,(p) = log(l + Vl - P 2 ) - Vl - P 2 + 1 - log2 = Pl . (19) 



It turns out that K{p c ) > for < p% < <fio(l) = 1 — log 2 = 0.3068... Essentially the pressure is given by the ideal gas 
contribution plus the first virial correction, with another correction which is exponentially small up to pi = 1 — log 2. 
Above the latter value a solution of Eq.s H15[) could not be found. 

These results are strongly consistent with the results of |3l| where it was shown that the resummation of the ring 
diagrams gives exactly Eq. I|17|) up to pi = 1 — log 2. At this value of the density a pole develops that seems to 
correspond to a liquid instability (the Kirkwood instability |2l|). 

It not clear if this instability implies that there are no solutions of Eq.s Ijl5(l (or of the HNC equations) for high 
densities, or it simply implies that one has to look for more complex solutions. Although the question is very interesting 
from the mathematical point of view (maybe also in relation to the problem of finding the most dense lattices 33] ), 
this question is not physically relevant in this contest. We will show later that the glass transition indeed preempts 
this instability that is therefore in a non-physical region of the density: this system becomes unstable toward replica 
symmetry breaking before reaching the Kirkwood instability. 

Using the exact relation 

P 1 dS 

- = l + -pY = -p— , (20) 



p 2 dp 



it turns out that 



\p for p < 1 

h 



S(p) = l-logp- I * , r ^„_™ B . 'I/ "V (21) 



\p(l + L{ Pc )e~ 2NK ^) for p > 1 , 



where, recalling that p^- = = N4> } [pa] , L(p) is such that 

L'{p)~2NK'{p)L{p) 1 1 

L[p)+ — — =1 * L{p) ~i-2K>( P w n (p) = 2 i + ^ -• (22) 

Again, up to exponentially small corrections, the entropy of the liquid is given by the ideal gas term plus the first 
virial correction. Similarly, by comparing Eq.s (|20|l and l|17|l we find that Y = 1 up to exponentially small corrections. 

V. GLASS TRANSITION IN LARGE SPACE DIMENSION 

To locate the glass transition we substitute Eq. 121|) in Eq. I|14|) and compute the equilibrium complexity 0, ITU 

HUSH!: 



E(p) = m 



2 Wm) 



dm 



= S(p) - I log[2^4*(l)] = S(p) - S mb (p) , (23) 
m=l z 
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FIG. 1: The function £(pi), Eq. (12 "il . for N = 100. Our computation is based on Eq. I|14|l which comes from an expansion in 
powers of the cage radius and is valid only close to the Kauzmann density where E = 0. The dashed part of the curve is then 
unphysical. Unfortunately, the exact point where Eq. 114H breaks down could not be estimated. Thus in the figure the point 
where the curve becomes dashed has been chosen arbitrarily. 



where S V ib(p) = | log[27rA*(l)] is the vibrational contribution to the entropy of the liquid; the glass transition density 
(or Kauzmann density) pk is defined by S(/o) = 0. Note that the cage radius is, from Eq. I|14fl . >/A oc p _1 ; thus for 
pi < it is exponentially large in N and the small cage expansion does not make sense. For pi > we get, neglecting 

the terms related to K(p c ) which are exponentially small in N, and using Vd ~ (^S-) 

1 



S(p) = 1 - log/9 - |p ~ 1 - Npx + N log(eTr) - N log N - log( V2ttAT 



N I nV2ttN' 
1 



3/2n 



,Npi 



/27T 



S vib (p) = d\og^—+ dlog[p(l + e - 2NK( ^)} ~ 2iVlog 



2tt 



(24) 



0.638 



0.638 



(27V + 3)7Vpi 



/27T 



and the equation for px is, neglecting for simplicity the terms growing slower than N, and defining a — 2 log q 638 
log(eTr) - 0.59, 

= £(/5i) Na- N log N + 2N Pl (N + 1) -^e Npi . (25) 

A plot of S(pi) is reported in Fig.^ Let us first neglect the terms proportional to N; then the equation S = gives 

1 



N log N = 2iV>i 



(26) 



(27) 



and a solution is simply Npi = ^ logiV. However the derivative of X evaluated in this solution is 

T,'( Pl ) = 2N 2 -^Ne Npl 

and is positive, thus this value corresponds to the unphysical zero of E (see Fig. QJ. The maximum of S is given by 
£' = 0, i.e. Npi = log4/V. Thus we look for a solution of Eq. (gSJ of the form 

e^ 1 = qN , (28) 

and we expect that q diverges with N. Substituting in Eq. I|26|l we get 

q = 41og(gV r /V) = 21ogiV + 41ogg , (29) 

which is solved iteratively by 

q = 2 log N + 4 log(2 log N + 4 log q) = 21ogiV + 41og[21og7V + 41og(21ogiV + 41ogg)] = ... . (30) 
Thus the value of the Kauzmann density is 

p K = 2N log N + 0(N log log N) , 



p lK = ilog(2AriogiV + 0(ArioglogiV)) 



(31) 
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A. Correction to the Kauzmann density 

For future convenience it is useful to compute the correction to ~p K due to the terms O(N) we discarded in Eq. 
The full Eq. is 

-N log N + 2N 2 Pl - ^e Npl =Na- 2N Pl . (32) 

We look again for a solution e Npi = Nq with q large. Then we obtain 

41ogg + 21ogA-g = 2a-— log(gA) . (33) 

The right hand side is o(l) while the left hand side diverges so the leading solution is q = qo with qo given by Eq. i|30|l . 
We look for a solution q = qo + qi with q\ <C qo : then the left hand side gives 

A\og(qVN) - q = 4 log (l + ?f) - qi ~ qi (--l) , (34) 

then 



9o/ V qo 



q 1 [ — -l)=2a + 0{N- 1 logN) 1 (35) 

qo 



and finally as qo ~ log A 
The result for ~p K is then 

with q given by Eq. (|30|) . 



-2a + 0((logiV)- 1 ) . (36) 



7>*r = ATgo-2Aa + o(-^ ) = Xij„ - . (87) 



B. Random close packing density 

F or p > p y the system is in the glass phase and the value to* that optimizes the free energy (per replica) is m* < 1 
0E 111 El to* is the solution of £(p,to) = 0, where 

d(<blm\ O' 
E(p,to) = to 2 V 1 J =S(p)-N\og(2n) + (2N + 3)N Pl -Nm-Nm(l-m)^f- 

+ A log to - 27V log(l - to) + 2 A log Q TO . 

As Qm ~ yJn/Am for m — ► (see Appendix^}, one can see from Eq. (|14|) that the cage radius A* — ► for to — * 0. 
The random close packing density p c is then the value of p at which to* = |19|. Then to compute p c we have to 
solve £(p, 0) = 0. Using Q TO ~ y^r/iro, 

EfaO)=S(p) + (2iV + 3)tfpi+ivQ-31og2j . (39) 



The condition E(p, 0) = 0, using Eq. (J23J, is similar to Eq. 



-A log A + 2A 2 pi - ^e^ 1 = JVa' - 2 Api , (40) 



with a' = 3 log 2 — \ — log(e7r) ~ —0.56. The solution is given by Eq. I|37|l with a — > a': 

p c = Ago - 2Aa' + O (r^j = ^qo + 1.12JV = p^ + 2.3A , 



Pic = i log(A ?Q + 1.12A) 



(41) 
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C. Lindemann ratio 



To check the consistency of the small cage expansion, it is interesting to estimate the Lindemann ratio in the glass 
phase, when p ~ 2N\ogN ~ dlogd. The Lindemann ratio L for a given solid phase is the ratio between the typical 
amplitude of vibrations around the equilibrium positions and the mean interparticle distance. In our framework it 
can be defined as 

L = P x l d ^fA (42) 
so that, using \f~A ~ 1/p from Eq. I|14|l. and p = pVd ~ d\ogd : one has 

y/dlogd ' ^ 
which is consistent with the assumption that vibrations are very small. 



D. Is the densest packing amorphous in large d? 



We can compare our prediction for the maximum density of amorphous packings p c ~ d log d with the best available 
bounds on the density of crystalline packings. The corresponding packing fraction scales as 4> c ~ 2~ d d\ogd. Unfor- 
tunately, the best lower bound for periodic packings is the Minkowski bound p ~ 1, while the best upper bound is 
the Blichfeldt's one, p ~ 2 d / 2 pfll Iscj. Our result for p c lies between these bounds so we cannot give an answer to 
the question whether the densest packings of hard spheres in large d are amorphous or crystalline. Hopefully better 
bounds on the density of crystalline packings will address this question in the future. However it is difficult to escape 
to the impression that the values of the densities of crystalline laminated lattices up to d=50 suggest that there 
are lattices where pi goes to a non-zero limit for infinite d. It is however quite possible that this is a preasymptotic 
effect. It would be very interesting to find the density of laminated lattices in larger dimensions. 

It would be also interesting to find out the maximum value of the density for which the inequalities (|15|l have a 
solution, allowing also delta functions in h(p): those equations are valid also for regular lattices; however the problem 
is hard and it is not clear how to attack the problem mathematically (it is a linear programming problem in an infinite 
dimensional space also for finite dimensions). 



VI. CONCLUSIONS 



Making use of the results of |l9l l3ll |32( we showed that the small cage expansion of the replicated free energy in 
large space dimension predicts a 1-step replica symmetry breaking transition. We obtained the asymptotic behavior 
of the glass transition density and of the random close packing density as d — > oo. It is worth to note that all the 
phenomena related to the glass transition happen is a region of densities which is very close to p = 1 , in the sense that 
pi = (d/2)^ 1 log pis bounded by o(logd/d), while the Kirkwood instability [3lJ,|32j happens for pi = 1— log 2 = 0.30 . . ., 
i.e. well beyond the interesting range of values of p\. In other words, the region of densities where the Kirkwood 
instability happens is never reached due to the glass transition happening at lower density. 

It would be interesting to estimate the corrections at finite d to the asymptotic expressions. For the simplified 
model (|15|l the corrections are exponentially small in d as found in . It is then reasonable that the exact expression 
for the liquid entropy as well differs from the asymptotic solution of l|15(l by exponentially small terms, and the same 
might happen for the replicated free energy. If this is the case, the results we presented in this paper should be exact 
for d — ► oo . We hope that future work will address this point. 
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APPENDIX A: THE FUNCTION Q(r) 
The function Q(r) has been defined in Eq. © as: 

Q(r) = J du 1 du 2 p(u 1 )p(u 2 )f(r + ui,u 2 ) - F(r) , 



(Al) 



where F(r) = -6(D - |r|) and f(r + Ui,u 2 ) = -1 + H a e~^ r+Ula ~ U2a \ Recalling that / dup(u) = 1 we can rewrite 



Q(r) = J du 1 du 2 p(u 1 )p(u 2 )l[e-^ r +^- u ^ - 8(\r\ - D) = F (r) m - 8(\r\ - D) 



where 



Fo(r) = / duidu 2 



(»ir-K-2) J 

e " 2A 



\r + u\ — u 2 \ — D)= I du 



e 4-4 



\r + u\-D) 



(A2) 



(A3) 



From the expressions above we expect that Q(r) is non zero only if \r\ ~ D. 
Let us first discuss the expansion of Q(r) in d = 1. Defining 



1, 



we have 



erf(t) = -= / dse -11 , 9(t) = -[1 + erf(t)] = -= / dxe~ x , 



fi (T — D 

du —^= 9{r + u-D) = 6 



(A4) 



(A5) 



Note that for small A, u is small too, and as \r\ ~ 13, the sign of r + u is the same as the sign of r. Thus we can write 
6(\r + u\-D)~e(r + u-D) + 6(-r -u-D) and 



F (r) = 9 



r - D 

VIA 



e 



' v/4l 



e 



V£3 



where in the last step we neglected contributions of order exp(— D 2 /A) for A — > 0, and 

M — Z3 N 



Q(r) = 



e 



f4A 



-D) 



(A6) 



(A7) 



From the expression above it is easy to see that Q(r) is non zero only if \r\ — D ~ yC4, as expected. This 
is because the function €)(t) is exponentially close to 9(t) if t is large. Given a smooth function one has 

/ drQ(r)f{r) oc \fAf{D) for small y/A, i.e. Q(r) cx \/A<J(|r| - D). We have now to compute the proportionality 
factor. Defining the reduced variable t = (r — D)/v4A, we have: 

(A8) 



(A9) 



drQ{r) = 2VA J dt[e{t) m - 6(t)] = 2VAQ{A) ~ 2 V AQ m + o{VA) 

~ V4.A 

where the function Q m is the limit for A — > of Q(A) and is given by 

/oc 
^[9(tr-e(*)] . 
-oo 

It is easy to show that Q m is a finite and smooth function of m for that 

Q m = (l-m)Q + O[(m-l) 2 ] , 

/OO 
dte(t) log 6(0 - 0.638 , 
-OO 



(A10) 



and that Q m diverges as Q m ~ y/n/4m for m — > 0. Thus we get the result Q(r) — 2^f~AQ rn 5(\r\ — D) in d = 1. 
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In dimension d > 1 we have, recalling that Fo(r) is rotationally invariant, for R = \r\: 

JdrQ(r)= JdrlF Q (r) m -9(\r\-D)}=n d J~dRR d - 1 lFo(Rr-d(R-D)}, (All) 
where f2d is the solid angle in d dimension, Sl^ = 27r d / 2 /r(d/2). The function Fq(R) can be written as 

FoiR) = I du \ 7 s^ m7+ul ~ Dh (A12) 

where i is the unit vector e.g. of the first direction in M. d . For small y/A, the u are small too, and the function 

6(\ri + u\ — D) is constant along the directions orthogonal to i. Thus we can show that the integral over the variables 
u M , fj, 7^ 1 is equal to 1 up to corrections of the order of exp(— D/\f~A). We finally get: 



f; i( /?!=/ du 1 ^=e(\R + Ul \-D) = e(^ 7 =^- ) , ( aui 



as in the one dimensional case. Again, the function Fo(R) m — 9(R — D) is large only for R — D ~ v A so at the lowest 
order we can replace with _D d_1 in Eq. IjAlljl . Finally, using that Fq(R) is given by Eq. jA"l3|) as in d = 1, we 

get 

J drQ(r) = (laD*- 1 J dR [F Q (R) m ~ 8(R - D)} = V d {D)2yfAQ m , (A14) 

where E<j(.D) is the surface of a d-dimensional sphere of radius D, T,d(D) = Q, d D 1 , i.e. in any dimension d we have 

Q(r) = 2y/AQ m d(\r\ - D) . (A15) 
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